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ABSTRACT 

Our analysis of Cassini/VIMS near-infrared spectra of Saturn’s Great Storm of 2010-2011 reveals 
a multi-component aerosol composition comprised primarily of ammonia ice, with a significant com¬ 
ponent of water ice. The most likely third component is ammonium hydrosulfide or some weakly 
absorbing material similar to what dominates visible clouds outside the storm region. Horizontally 
heterogeneous models favor ammonium hydrosulfide as the third component, while horizontally uni¬ 
form models favor the weak absorber. Both models rely on water ice absorption to compensate for 
residual spectral gradients produced by ammonia ice from 3.0 /mi to 3.1 /tin and need the third com¬ 
ponent to fill in the sharp ammonia ice absorption peak near 2.96 /mi. The best heterogeneous model 
has spatial coverage fractions of 55% ammonia ice, 22% water ice, and 23% ammonium hydrosulfide. 
The best homogeneous model has an optically thin layer of weakly absorbing particles above an opti¬ 
cally thick layer of water ice particles coated by ammonia ice. This is the first spectroscopic evidence 
of water ice in Saturn’s atmosphere, found near the level of Saturn’s visible cloud deck where it could 
only be delivered by powerful convection originating from ~200 km deeper in the atmosphere. 
Subject headings: Saturn; Saturn, Atmosphere; Atmospheres, composition; Atmospheres, dynamics 


1. INTRODUCTION 

Beginning in early December of 2010 amateur as¬ 
tronomers around the world observed the outbreak of 
a rem arkably l arge storm in the atmosphere of Saturn 
dSanchez-Lavega et all [Mil ). This large feature is the 
sixth known of a class of features historically referred to 
as Great White Spots, which ha ve occurred at ro ughly 
thirty year intervals (iSanchez-Lavega et al.111991D . i.e. 
at intervals close to one Saturnian year, always in the 
northern hemisphere, and usually in that hemisphere’s 
summer. This is the first such storm observed close up 
by a spacecraft (Cassini) in orbit around Saturn. The 
Cassini Imaging Science Subsystem (ISS) captured im¬ 
ages of the storm beginning in December (Fig. [l]) and 
the Visual and Infrared Mapping Spectrometer (VIMS) 
captured its first spectra of the storm in February 2011 
(Fig. [2]). It was then centered at 35° N and covered 
about 7° in latitude, with a long and widening down¬ 
stream plume. The vertical and horizontal scale of its 
effects were probed by thermal infrared spectroscopy, 
which showed that within a month the storm generated 
large perturbations of atmospheric temperatures, winds, 
and co mposition over a wid e area between 20° N and 
50° N dFletcher et alJl2 011bD. Dynam ical models of such 
storms (IHueso and Sanchez-Lavegal I2004D suggest that 
they originate at the 10-12 bar level where water vapor 
condensation is expected, and generate enormous con¬ 
vective towers reaching pressures as low as ^150 mb. It 
is thus conceivable that these storms provide a window 
into the deeper layers of Saturn’s atmosphere. 

The near-IR spectral observations of VIMS (from 1 
/tm to 5.15 /tm) are of special interest because they are 
particularly diagnostic of cloud particle composition. It 
was immediately obvious from VIMS spectra that the 
cloud particles in this storm had an unusual composition; 


the y disp l ayed a su rprisingly strong absorption near 3 /tm 
dBaines et al.ll20fTl) that was not seen anywhere else on 
Saturn (Fig. 2). Ammonia ice (NH3) and ammonium 
hydrosulfide (NH4SH), as well as water ice, are possible 
contributors because they are all strong absorbers in this 
region of the spectrum. While all of these molecules are 
also expected to produce clouds on Saturn, a thick over- 
lying haze has largely shielded them from view. Until re¬ 
cently, there has been little spectroscopic evidence of any 
of them except in association with small storm clouds in 
the southern hemisphere, where excess absorption near 
2.7 in n sugges ted the presence of freshly condensed NH 3 
dBaines et alj|2009l ). The Great Storm of 2010-2011 is 
by far the largest 3-/im absorbing feature ever seen on 
Saturn and provides the best opportunity for detailed 
spectroscopic analysis. 

Our objective here is to model the composition and 
vertical structure of the storm feature, using detailed 
spectral characteristics to determine the dominant com¬ 
ponents of the storm cloud particles, as well as to esti¬ 
mate the size of the particles and whether they are likely 
to be composite or pure. Of all the pure substances we 
tried, particles of NH 3 provided the best fit to the ab¬ 
sorption feature. However, pure NH 3 presents problems 
because its narrow strong absorption feature at 2.96 /tm, 
which should be resolved by VIMS, is seen only weakly 
if at all. Other characteristics of NH 3 do appear to fit 
well, but some mechanism that suppresses the 2.96-/tm 
feature seems to be required. Contributions of other ab¬ 
sorbers such as NH4SH and water ice can provide im¬ 
proved fits, suggesting that all three substances might be 
present in the cloud particles, with NH 3 being the most 
prominent contributor and water the most probable sec¬ 
ondary component. In the following we first present the 
VIMS observations, describe our spectral analysis tech¬ 
niques, evaluate spectral fits for pure substances, then 
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Fig. 1.— Cassini ISS images of Saturn’s Great Storm of 2010-2011 (A) at an early stage of development in December 2010 and (B) at 
approximately the same time as the first VIMS spectral images obtained in February 2011. Both images were made with a CB2 filter (752 
nm) and the wide angle camera. 
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Fig. 2.— (A-C) VIMS spectral images of the storm in February 2011, sampled at wavelengths indicated within each panel, remapped to 
an orthogonal projection, and linearly stretched so that white corresponds to a reflectivity (I/F) of 0.4. The color composite (D) is created 
from spectral images at 3.048 pm (blue component), 1.887 pm (green component), and 4.081 pm (red component). Numbered boxes 
indicate the locations at which spectral samples (E) were obtained. Here closely similar spectra are averaged to simplify the figure. Strong 
absorption at 3.048 pm reduces the blue contribution in the color composite image, making the storm feature appear orange. The parts of 
the storm that are bright at 4.081 pm (locations 1 and 2, for example) are much darker at 3.048 pm than the region ahead of the storm 
(locations 3 and 4). 1.887 pm and 4.081 pm are pseudo-continuum wavelengths, subject to little gas absorption, and thus show reflective 
properties of clouds more directly. The strong reflectivity at 4.081 pm indicates that the cloud particles are relatively large (~1 pm or 
more). The spectra from regions 3-6 are very similar, especially in the 2.9-3.04 pm region where they contain a similar phosphine absorption 
feature. The difference between those spectra and the storm head spectra (1 and 2) is colored yellow to show the extra absorption that is 
mainly due to particulates. 

consider alternatives involving combinations of NH 3 with as composite particles, and also as layers of different com¬ 
other substances, both as co-mingled pure particles and positions. 
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2. OBSERVATIONS. 

We used publicly available spectra acquired by VIMS 
during 24 February 2011 at 00:36:36 UT, data cube 
number 1677201862J3, acquired under observation name 
VIMS_145SA_WIND5HR001, at a phase angle of 52° and 
a range of 1.84xl0 6 km. Th e VIMS inst ru ment and 
investigation are d escribed by I Miller et all (|1996t l and 
I Brown et~a l. (2004). The instrument’s spectral range is 
0.35-5.1 /zm, with an effective pixel size of 0.5 milliradian 
on a side and a near-IR spectral resolution of approxi¬ 
mately 15 nm (sampled at intervals of approximately 16 
nm). The atmospheric pressures that can be sampled by 
this wide spectral range are indicated in Fig.[3l which was 
computed for a clear atmosphere and thus over-estimates 
the penetration depths below the 300 mb level, where sig¬ 
nificant aerosol opacity is a limiting factor. Methane is 
the most prominent absorbing gas at most wavelengths. 
Hydrogen collision-induced absorption is dominant in the 
1.95-2.15 /zm range. Phosphine is prominent in the 2.8-3 
(im and 4.2-5 /im regions. Arsine and CH 3 D have much 
smaller effects, and germane has too small an effect to 
be noticed in our investigation. 

A selection of these earliest VIMS observations of the 
Great Storm are shown in Fig. [2j Remapped images are 
shown at three key wavelengths and spectra from seven 
representative locations are displayed, including means 
of two just upstream of the storm head and two within 
the storm head itself. Compared to spectra at other lo¬ 
cations, the cloud reflectivities in the storm head are the 
highest at 4.081 /zin and the lowest at 3.048 //in. As 
shown in Fig. [2jy spectra obtained upstream (locations 
3 & 4) display absorption between 2.8 and 3 /zm that is 
due to phosphine, but show little absorption at 3.05 //m. 
The upstream clouds are also more transparent at ther¬ 
mal wavelengths, evident from their large apparent I/F 
values near 5 //m. The storm particles are brighter than 
the upstream cloud particles in most window regions, 
somewhat more absorbing at 2.7 //m , but dramatically 
more absorbing at 3.048 /zm. The image at 1.887 /zm 
(Fig. [2K) is used for the green component of the com¬ 
posite image because it provides the closest match to the 
I/F at 3.05 /zm for those regions that do not contain a 3- 
/zm particle absorption, both of which, along with 4.081 
/zm, are pseudo continuum wavelengths that are rela¬ 
tively unaffected by gas absorption. This means that the 
ratio of I/F at 3.048 /zm to that at 1.887 /zm should be a 
good measure of the degree of 3-/zm absorption, as illus¬ 
trated in Fig. [4] The regions of high reflectivity at 4.081 
/zm are places where the cloud particles are relatively 
large (probably 1 /zm or more), and strong absorption at 
3.048 /zm (Fig. 03) appears as bright orange in the color 
composite (Fig. HP). 

3. VIMS INSTRUMENT CHARACTERISTICS AND 
CORRECTIONS 

Before proceeding with the quantitative analysis of the 
VIMS spectra, we first describe special characteristics of 
the VIMS instrument and summarize how we handle data 
processing, photometry, and navigation. 

3.1. VIMS responsivity errors and corrections 

The responsivity of the VIMS detector array is signif¬ 
icantly reduced at joints between order sorting filters, 
which occur at wavelengths of 1.64 /zm, 2.98 /zm, and 


3.85 /zm (|Miller et al.ll 19961 1 Brown et al.ll2004ll . The de¬ 
cline is most severe at 1.64 /zm and 3.85 /zm, and VIMS 
spectra in t hese reg i ons ca nnot be trusted, as shown by 
ISromovskv and Frvl (j2010bl) for the 1.64 /zm joint, and 
by our own radiation transfer modeling for the 3.85 /zm 
joint. The spectra near the 2.98-/zm joint are reasonably 
well corrected in spectra obta ined during the Jupit er en¬ 
counter, as shown bv ISromovskv and Frvl (|2010bl) , but 
minor problems are present in the VIMS spectra of Sat¬ 
urn. The 2.98-/zm joint is unfortunately located so close 
to the prominent NH 3 ice feature near 2.96-/zm that its 
effects can be misinterpreted as being due to the pres¬ 
ence of NH 3 ice. In fact, all the raw VIMS spectra of 
Saturn do contain a dip very similar to what might be 
produced by a small component of NH 3 ice in the main 
tropospheric haze layer. However, the spatial variation 
of the feature is much more easily understood as an in¬ 
strumental artifact. The feature appears to be a frac¬ 
tional effect with no dependence on latitude or view an¬ 
gle, but only on the x-coordinate of the VIMS spectral 
data cube. It is not plausible that an atmospheric effect 
could produce this characteristic. Our interpretation is 
confirme d by ISO measu r ement s of Saturn’s 3-/zm spec¬ 
trum by lEncrenaz et alJ (|1999( 1 , which provide no evi¬ 
dence of such a feature. 

To correct for the responsivity error caused by the 2.98- 
/zm filter joint we collected spectra covering regions with¬ 
out strong convective features, then interpolated across 
the artifact for each spectrum and used the fractional 
depression below the interpolated line as a measure of 
the responsivity error that needed to be corrected. We 
found essentially no variation along image columns (y- 
direction), and thus used an average fractional depres¬ 
sion along each column to define the correction for any 
position along a column. We did find a significant vari¬ 
ation perpendicular to the column direction, however. 
The fractional error for band 127 varied smoothly from 
0.75 at column zero to 0.52 at column 63, a change that 
is about 30 times the standard deviation along that col¬ 
umn. A physical explanation for this behavior in terms 
of instrument design and operation remains to be deter¬ 
mined. With this correction we are able to separate the 
effects of the filter joint from those of true particulate 
absorption in the same region of the spectrum. 

3.2. VIMS errors at low signal level 

Althoug h a spec if ic cau se has not been identified, 
ISromovskv and Frvl ( 2010bl ) showed that very dark re¬ 
gions of VIMS spectra tend to be significantly brighter 
than independent observations from the ground and from 
HST using NICMOS. This occurs at very low signal lev¬ 
els, for which the instrument is responding with 1 Digi¬ 
tal Number or less. What appears to be a small offset, 
just a fraction of a DN, if not corrected for, leads to ex¬ 
cessive stratospheric haze opacity in models designed to 
match these spectra. As we have no well-defined correc¬ 
tion for this effect, it limits our ability to constrain the 
stratospheric haze, which is only noticed at wavelengths 
for which the atmosphere is extremely dark because of 
strong gas absorption. 

3.3. VIMS wavelength corrections 

The wavelength calibration tha t worked we ll for mod¬ 
eling VIMS Jupiter spectra (ISromovskv and Frvll2010bt f 
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Fig. 3. — Penetration depth of near-IR photons indicated by pressures at which a unit albedo reflecting layer produces an external I/F 
of 1/e at normal incidence and viewing. Left: the CCD spectral region. Curves are shown for methane only (gr een), ammon i a onl y 
(cyan), phosphine only (red), arsine o nly (blue), CH 3 D only (purp le), and all gases combi ned (black), assumi ng the Lin dal et al.l ((1985) 
te mperature profi l e, the He/H 2 ratio of Conrath and Gautier ( 2000), the ammonia profile of Prinn ct al] (| 19841) . and the phosphine profile 
of IFletcher et ahl (l2009al'l with a scale height of 0.375 times the gas scale height. Included, but not shown separately, is H 2 CIA, which 
peaks near 2.1 /mi. 



VIMS_145SA_WIND5HR001, 1677201862 

l/F ratio: 3.05 /xm / 1.89 /xm stretch=0-1.1 


quires shifts of the VIMS wavelengths that vary with 
wavelength. In Fig. [5J we show the results of match¬ 
ing 10 wavelength segments containing different absorp¬ 
tion bands. The required shifts are a significant frac¬ 
tion of the 15-nm VIMS resolution, ranging from 12 nm 
near 2 gm to -4 nm near 5 /im. We found that ramped 
wavelength shifts lead to better spectral fits than con¬ 
stant segments. The ramp functions we adopted, shown 
as dot-dashed lines, were adjusted to roughly follow the 
segment shifts , and t o cross zero near 5 /im, where the 
[Fletcher et al.l (|2011aD analysis of night-side Saturn spec¬ 
tra from 2006 found good agreement between VIMS and 
model wavelengths. As shown in Section 15.11 fit quality 
is much worse if unshifted wavelengths are used. This 
implies that the VIMS instrument changed behavior be¬ 
tween the early part of the 2000 Jupiter encounter and 
operations at Saturn in 2010-2011, by almost as much 
as the change from pre-launch to the time of the Jupiter 
approach. Whether the change in instrument tempera¬ 
tures between these two time periods could have been a 
factor remains to be determined. 


Fig. 4.— VIMS images of I/F ratio between 3.05 /im (numerator) 
and 1.89 /im (denominator). Absorption by NH 3 is greatest where 
the ratio is smallest. 

came from PDS example calibrated data cubes, which 
used a set of wavelengths that was different from those 
contained in the Jupiter flyby data cube headers. Ac¬ 
cording to iMcCord et al.l (12004T) . the VIMS wavelength 
calibration changed since launch by about 12 to 22 nm. 
Using comparisons with NIMS, they revised the VIMS 
wavelength calibration to what we believe is in the PDS 
example calibrated data cubes. These wavelengths are 
the same as appear in current PDS data cube headers for 
Satu rn (except for the last wa velength). They also match 
what lCruikshank et al.l (l2010h refer to as the RC15 wave¬ 
length calibration. But, while these wavelengths worked 
well for the December 2000 Jupiter flyby data, hen used 
with the 2010-2011 Saturn data cubes, they produce con¬ 
flicts between measurements and model calculations re¬ 
garding the positions of gaseous absorption bands. To 
align VIMS spectra with these absorption bands, which 
are accurately known from line-by-line calculations, re- 



Fig. 5.— VIMS wavelength shifts obtained from matching in¬ 
dividual spectral bands within 10 segments (solid stair-step with 
gray error bands), and adopted ramp shifts (dot-dashed lines) that 
provide improved matches to complete spectra. 

3.4. Image processing and navigation 

VIMS image cubes were processed using the pipeline 
processing code downlo aded fro m the Planetary Data 
System, described by IMcCord et al.l (120041) . To navi¬ 
gate the images, we extracted spacecraft and target body 
states, and instrument pointing information from the in- 
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dex.tab files contained in the PDS archive volumes. The 
sun to target distance was retrieved from Cassini NAIF 
SP kernel files, bypassing the VIMS pipeline code, which 
used a fixed distance (9.033 AU). These data allow ac¬ 
curate transformation between image coordinat es and 
planet coordinates using SPICE toolkit software (lActonl 
I1996T) . We processed the VIMS data cube up through the 
“specific energy” product stage using pipeline software, 
and then computed the I/F using the correct sun-target 
distance and the solar spectrum provided in the PDS 
volumes. The VIMS calibration uncertainty is at least 
as large as the near-IR. sol ar sp ectral irradiance uncer¬ 
tainty of ~3%-5% (|Colina et al.l 1 1 996h . It appears to be 
less than 10 % u ncerta in, based o n win dow-r egion com¬ 
parisons made bv lSromovskv and Frvl (|2010bH with inde¬ 
pendent spectra of Jupiter. 

3.5. VIMS noise characteristics. 

For most applications at moderate signal levels, the 
random noise leve l o f VIMS is very sm all, nominally 
less than 1 DN. iSromovskv and Frv l (j2010b ) tried to bet¬ 
ter characterize the noise by comparing spectra in cubes 
that imaged the same region of Jupiter with as small a 
time difference as possible. They found measurements in 
the same channel to be within about 1.2 DN RMS for 
a single measurement, even for signal levels up to 2000 
DN or more. However, at signals 1 00 tim es sma ller, as 
found near 2.3 /irn, ISromovskv and Frvl (l2010bh found 
a much larger fractional noise, which can be charac¬ 
terized as an I/F offset noise, which they crudely esti¬ 
mated as ~5xl0 -4 from the on-disk spectra. Other pos¬ 
sible sources of uncertainty, though not strictly random 
noise, are scattered light inside the spectrometer, wave¬ 
length errors, line-spread uncertainties, and wavelength- 
dependent absolute calibration errors. 

4. RADIATION TRANSFER CALCULATIONS 

4.1. Atmospheric structure and composition 

The assumed composition of the atmosphere as a func¬ 
tion of press ure is displayed in Fig. [ 6 ] We used the 
iLindal et alJ (1 981/ 1 temperature structure between 0.2 
mb and 1.3 bars and dry adiabatic extrapolation to ap¬ 
proximate the structure at deeper levels. iLindal et al.l 
(1198511 used a He/H 2 ratio of 0.06/0.94 in deriving their 
profile, and noted that the profile could be scaled to ac¬ 
count for a different He/I R rat io. We used t he revised 
ratio of 0.135±0.025 due to lConrath and Gautier! (|2000H . 
but did not scale the Lindal et al. thermal profile because 
that would result in disagreem e nt wit h the thermal emis¬ 
sion spectrum. The lConrath and Gautier! (12000H revision 
was motivated by a disagreement between the pre-Galileo 
He/H 2 ratio for Jupiter and the accurate in situ measure¬ 
ments by the Galileo Probe. The pre-Galileo result was 
based on a combined analysis of IRIS spectra and the 
Voyager-measured radio occultation refractivity profile. 
In an attempt to avoid whatever error had occurred in 
that Jupiter analysis, Conrath and Gautier based their 
Saturn revision entirely on Voyager IRI S observations . 
A lower value of 0.08 was given bv lFouchet et al. (2009') 
based on Cassini CI RS spectra and Cass ini radio occul- 
tations analyzed by iGautier et al.l (1200611 . Fortunately, 
these uncertainties in the He/H 2 ratio have little effect 
on our analysis. As long as the derived thermal profiles 
remain in agreement with the thermal emission spectra, 


which is true for all of these results, we found that our 
derived cloud model parameters are not very sensitive to 
the He/H 2 ratio. Test cases for ratios from 0.06 to 0.135 
show that the derived cloud parameters don’t change 
more than their uncertainty limits. After trying a va¬ 
riety of values for the CH 4 volume mixing ratio (VMR), 
we found the best fits were generally obtained with the 
IFletcher et al.l (l2009bll value (4.7±0.2)xl0~ 3 , which cor¬ 
responds to a C H/t/H? ratio of 5.3x 10 ~ 3 . For CH 3 D 
we also used the IFletcher et al.l (l2009bll VMR value of 
3x 10 -7 . The most important variable gas is PH 3 and its 
verti cal p rofile need s to be adjusted to fit VIMS spectra. 
IFletcher et al.l (l2009aT) has shown that PH 3 is enhanced 
over the equator and depleted over belts and within hot 
cyclonic polar vortices. We found that a reasonably ef¬ 
fective parameterization was a fixed deep mixing ratio 
and an adjustable scale height defining its rate of de¬ 
crease for pressures less than about 600 mb. Our best 
fits outside the storm region were obtained with a deep 
m ixing ratio of 6.4x 10~ 6 , which is the global mean value 
of lFletcher et al.l (I2009ajl . and a scale height of 37.5% of 
the pressure scale height. Inside the storm region bet¬ 
ter Gts were obtained with a more extended PH 3 profile 
using the same deep mixing ratio and 50% of the pres¬ 
sure scale height. The greater scale height within the 
storm is consistent with increased upward transport of 
PH 3 , as expected in regions of st ron g convectio n. We 
used the NH 3 profile given bv lPrinn et al. (1984), which 
has a constan t deep mixing ratio of 2.05 xl0~ 4 , based on 
iCourtin et alJ (119841) . More recent VLA results include 
a value of (4.8 ±l)xl0 -4 bv iBriggs and Sackettl (119891 1 
and some evidence for depletion of NH 3 in the 2-4 bar 
region, which they interpret as evidence for an NH 4 SH 
cloud and H 2 S vapor. Arsine has a noticeable affect on 
the VIMS spectra near 5 /.irn, which is where ammonia 
gas also plays a relatively minor role. Although NH 3 ab¬ 
sorbs over a wide wavelength range, its effects are largely 
shielded by overlying cloud opacity. It does however have 
a minor influence in blocking some of the thermal radia¬ 
tion at wavelengths beyond 4.6 pin. 


PH3 Model = F09sc.375 
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Fig. 6 .— Assumed volume mixing ratios of spectroscopically im¬ 
portant gases in the atmosphere of Saturn. See text for references. 
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4.2. Gas absorption models 

Because we wanted our correlated-k methane absorp¬ 
tion models to be based on line-by-line calculations, 
we limited our spectral analysis to wavelengths greater 
than 1.2 68 Atm. Recent adva nces in measuring line 
spectra dCampargue et al.1l2Q12H have improved line-by- 
line calculations so much that for A > 1.268 pm they 
are now likely more a ccurate than band models (see 
ISromovskv et ah ( 201 2 ;) for comparisons and correlated- 
k models). We used the same line data for comput¬ 
ing CH 3 D absorption models. For NH 3 we used the 
combin ed correlated-k abs orption model described by 
iSromovskv and Frvl f 2010aD . which is based pri marily on 
the Goody-Lorentz band model of iBowles et, al.l (|2008|1 . 
To model phosphine (PH 3 ) absorption we used our own 
ex pon ential sum appr oximations based on the line data 
of iButler et al.l (12006H in t he 2.8-3.1 pm region and on 
the iRothnran et al.l (l2009tl HITRAN 2008 line data in 
the 4.1-5.1 /xm region. Our model of A.SH 3 absorption is 
based on line data from lTarragol (1996) (via G. Bjoraker, 
via B. Bezard, personal communication). Where CH 4 
and PH 3 gas absorptions overlap we compute opacities 
for 100 combinations of 10 CH 4 terms by 10 PH 3 terms 
and then sort and refit to a 10 -term weighted sum, fol¬ 
lowing the suggestion of lLacis and Oinaa (|199lH . We do 
the same with all other gas overlaps in succession, which 
reduces all gas absorptions to ten-term exponential sum 
models at all wavelengths, even in regions with overlap 
of several gases. 

Collision-induced absorption (CIA) for H 2 and H 2 -He 
was calculated using programs downloaded from the At¬ 
mospheres Nod e of the Planetary Dat a System, which are 
documented by IBorvsowl (1991. 1993 ) for th e ll_> IF fun¬ 
damental band. IZheng a nd Borvsow (11995fl f or the first 
H 2 -H 2 overtone band, and bv IBorvsowl (119921 for H 2 -He 
bands. Although these programs and/or tables gener¬ 
ally provide coefficients only for normal and equilibrium 
hydrogen, absorption for any para fraction can be de- 
rived by appropriate li near combinations of these values 
(jBirnbaum et al.lll996H . We used equilibrium hydrogen 
generally and made test calculations with normal hydro¬ 
gen to assure ourselves that our conclusions were not 
sensitive to the ortho/para ratio. 

4.3. Multiple scattering methods 

For multiple scattering calculations we make use of a 
modif icat ion of the do ubling and adding code described 
bv ISromovskv! (j2005al fbl). The modification is to remove 
the Raman source term, which is not needed at near-IR 
wavelengths, and restore the or iginal b lack body emis¬ 
sion source term (lEvans and Stephens! Il99ll l to handle 
the 5-/tm region of the spectrum. We used a grid of 
44 pressure levels from 0.5 mb to 10 bars, distributed 
roughly in equal log increments, except that additional 
layers are introduced where cloud layers are inserted. 
As a function of wavelength we stop adding layers af¬ 
ter the one-way vertical gas optical depth exceeds six. 
To model the medium phase angle VIMS observations 
we used 16 quadrature points in both zeni t h angle and 
azimuth. We followed [Sromovskv and Frvl (j2010al lbl) in 
modeling coa ted-s phere particle s cattering using code 
originated by IT 0011 and Ackermanl (|198lH . We approx¬ 
imated the line-spread function of the VIMS instrument 
as a Gaussian of FWHM = 0.015 /xm, then collected all 


the opacity values within ±FWHM of the sample wave¬ 
length, weighted those according to the relative ampli¬ 
tude of the line-spread function, then sorted and refit 
to ten terms again. A special treatment is required in 
the 2.9-3.0 pm region where NH 3 cloud particles have 
very sharp absorption features. In this region model cal¬ 
culations were made at 5-cm _1 intervals (the maximum 
sampling frequency of our correlated-k models) and sub¬ 
sequently smoothed to VIMS resolution. 

4.4. Fitting cloud models to observations 

Our cloud models usually included both vertically ex¬ 
tended layers, where layer boundaries could be con¬ 
strained by the observations, and compact layers, where 
the vertical extent of cloud layers could not be con¬ 
strained. In some cases we evaluate both compact and 
vertically extended models, often finding that both can 
work for layers that are weakly constrained. For mod¬ 
els where the number of adjustable parameters is rela¬ 
tively small, parameters were adjusted to minimize y 2 
using a form of_the Levenberg-M arquardt algorithm, as 
described bv [Press et all (]1992f l. In most cases, model 
complexity is too large to reach automated convergence 
of all parameters in a practical time period (or ever). To 
reduce the scale of the fitting problem and increase al¬ 
gorithm stability, we use constraining assumptions for a 
subset of the parameters, which usually cannot be as¬ 
signed reliable uncertainty estimates. It is important 
to note that formal uncertainties in the remaining fitted 
parameters are only valid within the context of the con¬ 
straining assumptions and cannot be treated as absolute 
uncertainties. 

To speed the fitting process we reduce the total number 
of spectral comparisons to a set of 48 wavelengths that 
sample key spectral features and a wide range of gaseous 
absorption strengths. After fitting adjustable parameters 
to minimize y 2 for the sampled points, we then compute 
a complete spectrum for more detailed comparisons. To 
test for the presence of bias due to the specific wave¬ 
lengths sampled in constraining our fits, we did a fit to 
the Location 4 spectrum with the complete set of 225 
uniformly spaced wavelengths as constraints. The two 
resulting sets of fit parameters agreed within their un¬ 
certainties, indicating no significant bias. 

4.5. Estimating uncertainties in derived model 
parameters 

The constraining properties of the observations depend 
on the nature of measurement errors as well as on the 
errors in the modeling of atmospheric opacity. As previ¬ 
ously noted, the VIMS random I/F measurement errors 
in a given wavelength channel are negligible except at 
very low signal levels. Systematic errors are more of a 
concern, both in absolute calibration and in the wave¬ 
length scale, which can itself result in I/F errors when 
comparing to models at slightly incorrect wavelengths. 
A more important error source, which is not strictly ran¬ 
dom, but can vary from wavelength to wavelength, is the 
uncertainty in radiation transfer modeling due to the un¬ 
certainty in absorption coefficient models. The transmis¬ 
sion curves on which previous methane opacity model s 
were based had an assumed error of 0.05 (llrwin et all 
l2006f) . Although not well characterized, the correlated-k 
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models derived from line-by-line calculations also pro¬ 
duce uncertainties. 

To try to account for these effects we used the follow¬ 
ing procedure. Initial model fits are used to establish 
a rough characterization of the vertical opacity struc¬ 
ture of the atmosphere. We then calculate model spec¬ 
tra for the rough model, and for two perturbations: (1) 
optical depths offset by 0.01 and (2) optical depths in¬ 
creased by 10%. The resulting I/F differences are root 
sum squared with I/F errors due to wavelength uncer¬ 
tainties, and with an I/F offset uncertainty of 5xl0 -4 
and a measurement/relative calibration uncertainty as¬ 
sumed to be 1%. The I/F error associated with wave¬ 
length uncertainty is calculated from the derivative of 
the I/F spectrum with respect to wavelength times an 
estimated wavelength uncertainty of 0.002 /im, which is 
1/8 of a VIMS line width. In some cases we assume a 
simpler uncertainty model, which is not cloud model de¬ 
pendent, but is characterized by an I/F offset uncertainty 
of 0.002 and a fractional I/F uncertainty of 10%. The ex¬ 
pected y 2 value for the best fit is equal to the number of 
degrees of freedom Np (number of points of comparison 
minus the number of adjusted parameters) and the ex¬ 
pected uncertainty in the value is equal to y/2Np (a Icel 
Il995f) . However, both error models result in minimum 
y 2 values several times Nf. This means that either our 
error estimates are too low, or our physical models of 
the clouds are incomplete. There are many reasons for 
the latter to be true, including not-accounted-for effects 
of non-spherical particles, multi-coated particles (more 
than two compositional layers), and heterogeneous struc¬ 
tures, as well as more complex vertical distributions. To 
approximately correct for bad error estimates and in¬ 
complete physics, we rescaled our y 2 values by the fac¬ 
tor Nf/xm IN before finding confidence limits. Because 
higher than expected y 2 values are probably not simply 
due to bad error estimates, this procedure is not entirely 
valid, but is clearly better than making no adjustment. 

5. MODELING CLOUD STRUCTURE OUTSIDE THE 
STORM REGIONS. 

5.1. Sample fits using a simplified cloud structure 

Near IR spectra outside the storm regions can be 
very well fit by relatively simple cloud structures, which 
are similar to those used successfully t o fit spectra 
in the CCD range from 0.3 to 1 um (IMunoz et ah 
2004: iKarkoschka and Tomaskdl2005t iPerez-Hovos et al. 

20051) . These structures contain a vertically diffuse 
stratospheric haze, an optically thick and physically thick 
main cloud layer, which requires visible and UV absorp¬ 
tion, but can fit near IR spectra with conservative cloud 
particles, and a deeper layer, to which most wavelengths 
have little sensitivity, but which can be constrained by 
the thermal emission that can be sensed in the 5-/tm 
region. Fits of such a model structure to the VIMS spec¬ 
tra from all non-storm locations identified in Figs. [2] are 
summarized in Tabic [l] where parameter values, uncer¬ 
tainties, and x 2 values are all derived from the subset 
of 48 spectral samples discussed previously. In all these 
fits we characterized the deep cloud as physically com¬ 
pact, but optically thick (arbitrarily set to r = 30) and 
absorbing (set to w=0.9 to provide adequate blocking of 
thermal emission). The only adjustable parameter of this 


layer was its effective pressure. As will be shown in the 
following, other deep cloud options could also be made to 
work. Cloud particles in the main haze layer (generally 
between 120 and 700 mb) were treated as Mie particles 
with refractive index n=1.4+0i. We fit the upper and 
lower pressures of this layer and its optical depth at 2 
/tm, and within the boundaries of this layer we assumed 
a vertical scale height equal to the gas scale height. Be¬ 
cause our x 2 values are several times the expected value 
(the number of degrees of freedom), it is clear that ei¬ 
ther our model lacks needed physics or our measurement 
error estimate is too small, by a factor of \JXmin/^f, 
or some of both. If we assume that our error estimate 
is the problem, then we need to multiply our standard 
parameter uncertainty estimate (the change needed to 

produce Ay 2 = 1) by the same factor of i /xmin/^f- 
Error estimates in Table |T] have this correction applied. 

The complete best-fit model spectrum for location 4 
is displayed in Fig. [7] We also evaluated the fit qual¬ 
ity of this spectrum using y 2 computed for all measured 
points, assuming a simplistic random error model with 
a reflectivity (I/F) fractional uncertainty of 0.1 and an 
offset uncertainty of 0.002, which provides more of a rela¬ 
tive than absolute judgment of fit quality. This provides 
a convenient way to compare full resolution versions of 
the best-fit model spectra. Two y 2 values are shown 
in the legend of Fig. [TJ the one labeled y 2 /7V(TOT) 
is computed for 1.268-5.15 (im wavelength range, and 
the one labeled y 2 /V(ABS) is computed for the 2.65-3.2 
/tm range where particulate absorption features can be 
prominent (in the storm region). All are computed with 
respect to the wavelength-shifted spectrum. If computed 
relative to an unshifted spectrum the y 2 /N (ABS) in¬ 
creases from 0.26 to 2.3 and y 2 /N (TOT) increases from 
0.98 to 3.8. This is a dramatic change that clearly con¬ 
firms the need for wavelength corrections. Although the 
fit quality for this model is very high, the model is not 
unique. A deep cloud that is vertically extended, with a 
constant optical depth per unit pressure, can fit just as 
well, provided it contains some absorption that increases 
with wavelength. 

5.2. Sensitivity of model spectrum to model parameters 

Figure[7]also displays the effect of removing each model 
layer one at a time, and also the effect of removing 
various gases from the atmospheric absorption model. 
Next to methane, phosphine has the most significant ef¬ 
fect on the spectrum. Hydrogen CIA is next most sig¬ 
nificant, but limited to wavelengths from 2 to 2.2 pm. 
The stratospheric haze provides the baseline I/F in the 
darkest parts of the spectrum. However, because of the 
anomalously large I/F val ues measured by VIM S at very 
small I/F levels (iSromovskv and Frvl l2010bt l. and the 
small reflectivity of the small stratospheric haze parti¬ 
cles at near IR wavelengths, what can be learned about 
that haze from the VIMS observations is rather lim¬ 
ited. We found that a Mie particle radius of 0.1 /mi 
or a Henyey-Greenstein scatterer with wavelength depen¬ 
dence of A -3 5 , a pressure level of 2 mb, and an adjustable 
optical depth provides a suitable model. But this is not 
a very realistic physical model of the stratospheric haze. 
Extending the haze downward as prior CCD-based mod¬ 
els have done, does not improve the accuracy of fits here, 
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probably due to the offset errors of the VIMS instrument. 

5.3. Alternate vertical distributions 

The variety of vertical distributions that can fit the 
near IR spectra is similar to the variety that have been 
used to fit the visible spectra, as illustrated in Fig. [HI 
where we compare our fits to the n on-storm spectrum 
from l ocati on 4 in Fig. [2] to models of Perez-Hovos et all 
(|2005f) and IKarkoschka and Tomasko (1200511 applicable 
to similar latitudes, but different time periods. The dis¬ 
agreements are not surprising, given the secular varia¬ 
tion that might have occurred between the relevant ob¬ 
serving times. We tried both broadly diffuse vertic al dis¬ 
tributions suggested by IKarkoschka and Tomaskol (120051) 
a nd th e more compact models used bv lPerez-Hovos et all 
(2005'!). We found reasonably good fits were possible with 
both types of models, even to the 5-/^m region. By mak¬ 
ing use of thermal wavelengths (A > 4.8 n m), we were 
able to place better constraints on the deep cloud struc¬ 
ture than were the previously cited authors, though here 
again there is not a clear distinction that can be made 
between compact and diffuse models. If we assume that 
the upper cloud has a real index of 1.4, and that the 
deep cloud is nearly conservative (-07=0.98) and compact 
(Model B in Fig. ® , then we find that the cloud needs 
to reside near 3.3 bars and have a high optical depth 
(r ~ 150). But if the cloud is more strongly absorbing 
(o7 =0.9 as in Model C in Fig.®, then the optical depth 
requirement is greatly reduced (r ~ 30) while the pres¬ 
sure remains essentially unchanged. In all these cases 
the main haze layer visible from above was assumed to 
have a refractive real index of 1.4, to be uniformly mixed 
with gas, and to have adjustable top and bottom pres¬ 
sures, which reached typical values of 178 mb and 780 
mb respectively, while the optical depth did not deviate 
much from 6.8 at 2 pm. Although the compact lower 
cloud models with a separated upper tropospheric layer 
(extending from 178 mb to 780 mb) provided the best 
fits (B and C in Fig. ®, we were also able to get rea¬ 
sonable fits with a semi-infinite diffuse cloud extending 
downward from 178 mb, using an opacity of 11 optical 
depths/bar (A in Fig.®, with a real index of 1.4 and an 
imaginary index linearly increasing from 0.001 at 1 pm 
to 0.003 at 5.5 pm. This la tter model is similar to that 
used by IKarkoschka and Tomask<3 (120051) . Perhaps the 
most significant deviation is that our VIMS-based fits 
for the stratospheric haze layer are not consistent with 
that haze extending down to 100-300 mb, as found by 
prior fits. We believe that this disagreement arises from 
an artifact of the VIMS instrument in which very low 
I/F values are too high, as explained in Section [3~2l 

5.4. Latitudinal variations in non-storm regions 

Fits to spectra from all non-storm locations identified 
in Fig. ® and summarized in Table ® are graphically 
displayed in Fig. ® In these non-storm fits, which sam¬ 
pled latitudes from 20° N to 48° N, there is a trend for 
the main upper tropospheric haze layer towards smaller 
particles and smaller optical depths with increased lati¬ 
tude. In these fits there is no significant variation in the 
effective height of the thermal blocking layer. The sig¬ 
nificantly worse fit to the location 6 spectrum, which is 
from the near equatorial region, arises mostly from the 


low model I/F in the 1.5-1 .6 pm continuum region, which 
might be due to a small amount of unaccounted-for ab¬ 
sorption in the main cloud later in that broad equatorial 
band. Note that these fits exclude the latitude region 
between 22° N and 30° N, where there has been a strong 
storm influence, as indicated by the temporal changes ev¬ 
ident in Fig.® Analyzing the character of that disturbed 
region is left for future work. 

5.5. Composition of the main cloud layer in non-storm 

regions 

None of the fits to the non-stornr regions required any 
unusual absorption in the vicinity of 3 pm, where NH 3 , 
N 2 H 4 , and NH 4 SH have substantial absorption. Accord¬ 
ing to Fig. ® there is a wide spread lack of 3-pm par¬ 
ticulate absorption in regions not affected by the storm. 
It thus appears that the main upper cloud layers outside 
the storm-affected regions are not composed primarily of 
these materials. This is a well known mystery that re¬ 
mains to be definitively solved. It is quite amazing that 
the main cloud layer on Saturn that is accessible to re¬ 
mote observations, which is of high optical depth and 
thus represents a significant cloud mass, is apparently 
made of materials that have no detectable near-IR spec¬ 
tral signatures and thus do not seem to contain any sig¬ 
nificant contributions from the main condensable gases 
on Saturn. 

Explanations for how ammonia could really be the 
dominant condensable in the main haze layer, even 
though it has not been detectable outside of storm re¬ 
gions, have focused on three possible concealment mech¬ 
anisms: ( 1 ) large particle sizes, ( 2 ) non-spherical particle 
sh ape effec t s, or (3) foreign coating or contamination. 
IWest et all (12009D reviewed these explanations and re¬ 
jected them all, noting that a foreign contaminant would 
have to be the dominant cloud mass to be effective at 
suppressing the NH 3 signatures. Our calculations with 
composite particles are consistent with the conclusions 
of West et al. and our inferred particle sizes are also too 
small (~l/rm or less) to support ( 1 ). 

The leading candi date for the upper tropospheric haze 
according to lFouchet et al . (200 11) is diphosphine (P 2 H 4 ), 
which does not absorb much at 3 /irn, but has a distinc¬ 
tive double absorption p eak si gnature near 4.3 /iin and 
another at 6 /im (lNixonlll956f ). If the 4.3 -/im. peaks of 
diphosphine have roughly the same strength as the dou¬ 
ble absorption peak of hydrazine (N 2 H 4 ), then that sig¬ 
nature should be visible even though it is partly obscured 
by PH 3 absorption. However, though no such signature 
is seen either in or out of storm regions, it is premature 
to completely rule out diphosphine until refractive index 
measurements allow accurate radiative transfer calcula¬ 
tions to be made. We can rule out hydrazine as the 
major component of the main upper tropospheric haze 
because its distinctive features would be easily detected, 
as apparent from the discussion in the following sections. 

6 . MODELING THE 3-/*m ABSORPTION IN STORM CLOUD 
SPECTRA. 

We first consider what compositions might be plausi¬ 
ble, try simple spectral fits using pure substances in a 
single main cloud layer, evaluate their relative advan¬ 
tages, then consider more complex structures including 
overlying sub-layers and composite particles. 
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Fig. 7. — Example model spectrum (thin solid curve) fit to the VIMS mid-latitude spectrum (thick gray curve) from location 4 in Fig. [2] 
The model parameters are given in Table [I] Contributions from each layer can be inferred from the spectral changes produced bv removing 
it. Relative fit quality for the entire spectrum (excluding the two regions indicated by vertical gray bars) is indicated by the x 2 /N (TOT) 
values, and for the 2.65 pm - 3.2 pm region of particulate absorption by x 2 /N (ABS) values, where there are 33 points in the ABS region 
and 223 in total. Sensitivities to gas absorptions by NH 3 (dotted), PH 3 (solid), ASH 3 (double dot-dash) and normal vs. equilibrium H 2 
(dot-dash) are indicated by ratio spectra in the bottom panel (modified model/best-fit model). Greater sensitivity to these gases would be 
obtained in regions of lower cloud opacity. The importance of thermal emission in constraining deep cloud opacity can be seen from the 
large I/F increase caused by removal of layer 3. Vertical gray bars indicate where order sorting filter joints cause uncorrected VIMS I/F 
uncertainties. At the joint near 2.96 pm, the measured spectrum has been corrected for responsivity errors following a procedure described 
in the text. The thin lighter gray curve shows the original uncorrected spectrum. The good agreement of the corrected spectrum with 
model calculations in this region attest to the efficacy of the correction. 
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Fig. 8 .— Sample cloud structure models for Saturn according to lPerez-Hovos et al.l <120051) and IKarkoschka and Tomaskol <120051) and, to 
the right of the vertical dashed line, our model fits to the spectrum from location 4 in Fig. [ 2| Structure A uses a u ses a uniformly mixed 
single particle type from cloud top to the bottom of the atmosphere, as used bv IKarkoschka and Tomaskol <120051) . To provide sufficient 
absorption of lo ng w avelengt hs this requ ires an imaginary index that increases with wavelength. Structures B and C, which are similar 
to that used bv IPerez-Hovos et all (120051) , differ in the single-scattering albedo assumed for the lower cloud layer. Higher albedo requires 
larger optical depth to block thermal emissions, x 2 values shown are for fits to the 48-wavelength sampled spectra. 


6 .1. Optical properties of candidate cloud materials 

Cloud materials expected in Saturn’s atmosphere that 
also absorb light in the 3-//m region, include NH 3 ice, 
NH 4 SH (ammonium hydrosulfide), water ice, N 2 H 4 (hy¬ 


drazine), and possibly solid P 2 H 4 (diphosphine), al¬ 
though it appears th at the latter h as only very weak 
bands in this region (|Frankisslll968[ ). Refractive index 
plots for these materials (except diphosphine, which has 
not been characterized) are shown in Fig. 1101 Among 
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TABLE 1 

Cloud model parameter fits for February 2011 non-storm spectra. 


Parameter 

Location 3 

Location 4 

Location 5 

Location 6 

Planetocentric Lat., ° 

36.78 

34.41 

47.52 

19.00 

Planetographic Lat., ° 

42.57 

40.10 

53.32 

22.94 

East Long ., 0 

-116.09 

-115.89 

-110.63 

-102.90 

Stratospheric haze r x 10 3 

8.6T0.9 

8 . 6 T 0.8 

7.5T0.7 

15.Oil.8 

Main haze Ptop, mb 

143±11 

178T12 

132ill 

lllil 2 

Main haze Pbottom, mb 

844±110 

785i 82 

577± 76 

633i 73 

Main haze r 

6.53T0.81 

6.75T0.66 

3.48T0.36 

8.57i0.83 

Main haze part, radius, pm 

0.55i0.04 

0.61T0.05 

0.43T0.02 

0.71i0.07 

Deep cloud Pbase, bar 

2.98T0.11 

3.25T0.07 

2.68T0.15 

2.99i0.11 

X 2 

204.08 

149.87 

169.30 

430.19 


Uncertainties are 1 -a values after adjustment for excess x 2 minimum, as described in text, x 2 
values are for fits to the 48-wavelength sampled spectra. Locations refer to Fig. [2] 


Fit parameters outside of storm region for February 2011 
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Fig. 9.— Pressure boundaries of main cloud layer and effective 
pressure of thermal barrier cloud vs latitude, with annotated values 
of location, optical depth, and particle size. Locations refer to Fig. 
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the more plausible storm particle components, NH 3 ice 
has the strongest and sharpest absorption at 2 /irn (the 
2.25 //m absorption feature would not be visible due to 
overlying gas absorption). Hydrazine, which is not ex¬ 
pected to be very abundant on Saturn, has a pair of 
strong absorption peaks between 3 and 3.2 /im, which 
would be very apparent unless cloud particles were very 
large. Though diphos phine, which is expec ted to be more 
abundant on Saturn (IFouchet et al.lf2009h . could plausi¬ 
bly also contribute to the storm particle mix, we have 
no evidence in the VIMS storm spectra of its distinctive 
double absorption peaks near 4.35 /im, and no refractive 
index measurements that would allow us to include it in 
quantitative radiative transfer models. NH 4 SH has an 
absorption that is roughly comparable to that of NH 3 at 
3 /im, but lacks a sharp spectral feature and continues 
to increase beyond 3.1 /im while NH 3 absorption drops 
significantly. 

There is some disagreement about the detailed re¬ 


fractive index variation w ith wavelength fo r NH3 . A 
com parison of results fr om iMartonchik et all ( 198411 and 
from IRoux et all (119791) is provided in Fig. [TTJ Because 
iRoux et al.l ( 19791) provides better spectral resolution in 




Fig. 10.— Imaginary (A) and real (B) components of the refrac¬ 
tive index vs. waveleng t h for candidate 3-pm a bsorbers, including 
HenO results of l Warre nl ( 1984 ). NH 3 results of IRoux et alj (119790 
and IMartonchik et al.l ( 1984 ) , NH 4 SH results of Howett et al'l 
(120071) , and N 2 H 4 results of I Clapp and Miller! (119961) . Vertical 
dotted lines at 2.96 and 3.05 pm indicate approximate locations 
of NH 3 absorption features. 

the 2.7-3.1 fi m region and pro vides no informat i on be- 
low 2.8 /im, we initiall y used IMartonchik et al.l (119841 ) 
for A < 2.8 /im and IRoux et al.l (1 19791) for longer wave¬ 
lengths. These two references also disagree on the level 
of absorption near 5 /im. We later us e d a di fferent blend 
in which we used Martonc hik et al.l (I1984H everywhere 
except between 2.88 /tm and 3.14 /tm, where IRoux et al.l 
dl979h is used. This provides improved fit quality for 
some models. There is also evidence for small var iation s 
of optical properties with temperature dSill et al.lll980lh 
which we cannot account for and thus are forced to ig¬ 
nore. 


6 . 2 . Initial spectral fits with pure particles 

The five substances we considered in preliminary cal¬ 
culations with single-composition particles are these: ( 1 ) 
conservative substance with refractive index n=1.4+0i, 
(2) NH 3 , (3) NH 4 SH, (4) N 2 H 4 , and (5) water ice. As 
indicated in Fig. [10] substances 2-5 all have significant 
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Fig. 11.— Comparison of IMartonchik et al.l 11984 ') and I Roux et al. ; (1979) measurements of the refractive index of NH 3 . Our secon d 
blende d version (dashed curve), uses the former everywhere except bet ween 2.88 pm and 3 .14 pm , where the latter is used. IRoux et al.l 
lil979ft provides a better match to VIMS observations near 2.9 pm, while IMartonchik et all ( 1984 ) provides a better match near 5 pm. 


absorption in the 3-pm region. The question we are try¬ 
ing to answer is which substance is dominant, if any, 
and is there evidence of contributions from more than 
one substance? The great depth of vertical convection 
suggest ed by numerica l mode ls of the dynamics of this 
feature dHueso and Sanchez-Lavegall2004i ) makes it plau¬ 
sible to consider that even water ice might reach levels of 
the main cloud layer. Intense lightning strikes inferred 
from static electric discharge (SED) events measured by 
the Cassini Radio and Plasma Wave instrument (RPWS) 
occurred mainly near the head of the storm and indicate 
that this is a deep convective storm area energized by 
latent heat of water condensation in the 10-20 bar level. 
Other lightning strikes seen optically by the ISS instru¬ 
ment confirm their ve rtical location 125-2 50 km below 
the upper cloud layer (IDvudina et al.ll20Tol) . 

To gain a better understanding of how various plau¬ 
sible particle constituents might contribute to the 3-/im 
absorption feature, we first try to fit the storm spectra us¬ 
ing pure substances in a single middle cloud layer, leaving 
for later the issue of possible composite particles and/or 
multiple layers of different composition. We assumed 
that the middle layer was composed of pure spherical 
particles, then adjusted the optical depth, particle size, 
and top and bottom pressures to optimize fits for the re¬ 
gion from 1.268 fim to 5.15 /irn, but ignoring the region 
from 2.65 (im to 3.2 /im, where particulate absorption 
is most prominent, as well as ignoring the two narrow 
regions indicated in Fig. [T2l where VIMS measurements 
have large uncertainties (described in Section I3TT1) . That 
allows the vertical structure of the aerosols to be primar¬ 
ily constrained by methane gas absorption and not much 
affected by particulate absorption. We then compared 
the model spectra with measured spectra in the 2.65 fim 
- 3.2 /im region to determine which provided the best 


compositional match. 

The fits for the five pure substances are displayed in 
Fig. I12I V. The fit parameters and fit quality are summa¬ 
rized in Tables [2] and [U We evaluate the overall fit using 
X 2 computed as described for the non-storm regions, he 
conservative particles with n=1.4 (Fig. [T714) provide a 
reasonably good spectral fit at all wavelengths except 
in the 3-fim region, where the model has a small atmo¬ 
spheric absorption due to phosphine and the measure¬ 
ments require a much larger absorption over a broader 
range that must be due to cloud particles. The phos¬ 
phine absorption is particularly small here because the 
cloud layer containing these particles is relatively high 
(140-550 mb) and optically thick (r=9 at 2 /im). The 
residual absorption (colored yellow in Fig. [TTR). which is 
the difference between the conservative model (red curve) 
and the measured VIMS spectrum (thick gray curve in 
Fig. 1121) . thus must be provided predominantly by parti¬ 
cle absorption, and entirely so between 3.05 and 3.1 /im. 
This also shows that the particle absorption is significant. 

The ammonia spectrum provides by far the best match 
to the measured spectrum of the storm head, yielding a 
X 2 value in the absorbing region that is 3.5 times better 
(smaller) than that of the next best fit. While ammo¬ 
nia particles provide the best fit among these simplistic 
single-component models, the resulting model spectrum 
still contains substantial differences from the VIMS mea¬ 
sured spectrum, especially the sharp ammonia ice ab¬ 
sorption feature near 2.96 fim, but also the excess bright¬ 
ness near 3.1 fim. The model spectrum also has a small 
dip at 2 /im, where NH 3 has another sharp, but much 
weaker, absorption feature, while a corresponding dip is 
not seen in the measured spectrum. 

Turning to the option of pure NH4SH particles (green 
curve in Fig. fl2l4). we see that the overall fit is better 
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2 3 4 5 

Wavelength, nm 



Wavelength, \im 

Fig. 12.— (A) A storm head spectrum from (location 1 in Fig. [2j, and fits to that spectrum using a main cloud layer consisting of 
either conservative particles (red), pure NH 4 SH particles (green), pure NH 3 particles (blue), or pure H 2 O particles (cyan). The vertical 
gray bars indicate regions excluded from our analysis because order sorting filter joints lead to uncertain responsivity corrections at those 
wavelengths. Relative fit quality for each entire pure spectrum (excluding the two regions indicated by vertical gray bars) is indicated by 
the x 2 /N (TOT) values, and for the 2.65 fim - 3.2 ^m region of particulate absorption by x 2 /N (ABS) values, where there are 33 points 
in the ABS region and 223 in total. The specific fit parameters for each spectrum are given in Table [ 2 ] The yellow region between the 
conservative particle spectrum (red) and the VIMS measured spectrum is due to cloud particle absorption. (B) As in A except that the 
blue model spectrum is a linear combination of pure spectra , with fractions adjusted to minimize x 2 i n the particle absorbing (ABS) region, 
and the red spectrum is for Model 3 (discussed in Section I6.4.3|> , which is a multi-layer horizontally homogeneous model with the main 
upper tropospheric haze layer containing an upper sub layer of conservative particles and a lower sub layer of ammonia-coated water ice 
particles, as defined in Table [ 6 ] 


than that obtained with conservative particles, mainly 
due to the better fit in the 3-/xm region, but the fit qual- 
ity in the 2.65-3.2 /rm region (y 2 /N (ABS) = 88.4) is 
ten times worse than obtained from pure NH 3 . It does 
have a potentially useful characteristic, however. Its I/F 
at the edge of that window («3.15 /im ) is lower than 
obtained from NH 3 , so that combining the two might 
improve the overall fit quality. While hydrazine (N 2 H 4 ) 
provides a decent overall spectral fit, its double absorp¬ 
tion peaks in the 3-3.2 /im region are a poor match to 
the observed absorption and also provide little help as a 
minor contributor. 

The fit for pure water ice (cyan curve in Fig. |T2jA) 
shows that water is a very strong absorber, with a spec¬ 
tral signature completely different from that of the other 
substances, and especially different from what is needed 
to match the observed spectrum. However, it does turn 
out to have useful properties when we consider composite 
fits. We often find that good fits using NH 3 and NH 4 SH 
result in the window regions outside of the 2.65-3 /im 
window being relatively too dark. By adding a small 
amount of the water spectrum, we can brighten these 


window regions without brightening the 3-/tm region be¬ 
cause of water’s strong absorption at those wavelengths. 
Incidentally, it may appear strange that water ice fits 
the observed spectrum better than NH 4 SH according to 
y 2 values, yet visually appears to be a much worse fit. 
The visual appearance is misleading because the plots 
use a logarithmic y axis, while y 2 calculations use a lin¬ 
ear weighting of differences (squared). 

6.3. Linear combination spectral fits 

We next consider linear combinations of pure spec¬ 
tra to see if mixed particle compositions could provide 
a better match than pure particles. These combinations 
can be thought of as spatially heterogeneous mixtures, 
perhaps as different particles visible in spaces between 
clouds of other particles. They are also suggestive of 
what might happen if an optically thin layer of one cloud 
type would overlie a layer of another cloud type (but not 
vertically separated by large distances), or if a cloud con¬ 
sisted of composite particles, perhaps with a core of one 
(or more) materials with a condensed shell of another ma¬ 
terial. The best linear combination, shown by the blue 
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TABLE 2 

Adjustable fit parameters for different main cloud compositions. 


Main cloud 
composition 

Main cloud layer parameter values 

Top Bottom Particle 2-/xm Optical 

P, mb P, mb radius, /im depth (2-/im) 

x 2 

n=1.4+0i 

129±3 

647T40 

1.34±0.07 

12.7T0.7 

538 

nh 3 

133±3 

849±90 

1.35T0.06 

18.7T2.1 

426 

nh 4 sh 

136±3 

945T120 

0.61T0.01 

7.5T0.7 

591 

h 2 o 

182±5 

401±60 

0.83T0.10 

16.4±3.7 

1274 

N 2 H 4 

138±3 

826±84 

1.13T0.04 

13.8T1.5 

479 


Fixed parameters included the pressure (2 mb), and 2-^rni optical depth (0.01) of the 
stratospheric haze (assumed to be a conservative, isotropic scatterer with A -3-5 wave¬ 
length dependence), and the pressure (1.7 bars), optical depth (100), single-scattering 
albedo (0.9), and asymmetry parameter (0.6) of the deep cloud. Mie scattering was 
assumed for the main cloud layer, with equal gas and particle scale heights. Uncer¬ 
tainty limits given here are formal errors useful for relative comparisons only, with no 
correction to account for excess % 2 minima. values listed here are for fits to the 
48-wavelength sampled spectra. 


TABLE 3 

Fit quality for single-particle spectra and best fit combination spectra. 


n=l. 4+0i 

Component fractions 
NH 3 NH 4 SH 

h 2 o 

Fit quality (x 2 /lV) 
1.26-5.16 Atm 2.65-3.1 Atm 

Notes 

1 




46.3 

292.2 

1 


1 



3.3 

8.8 

1 



1 


16.0 

88.4 

1 




1 

16.3 

72.4 

1 


0.55T0.04 

0.23T0.02 

0 . 22 ± 0.02 

2.0 

1.1 

2 

0.13T0.02 

0.62±0.03 


0.35T0.03 

2.4 

2.9 

2 

O.OOiO.OO 

0.83T0.02 

0.17T0.02 


3.7 

5.2 

2 


1 Fits for single particle types in the main cloud layer, with parameters given in Table [2] 
2 These are linear combination fits with blank columns indicating particle types that were 
not included in the fit. Here x 2 values were computed using the simplified error model 
described in the text. Note the significant degradation in fit quality when water ice is 
excluded. 


curve in Fig. [T2t3. uses a mix of 55±4% NH 3 , 22±2% 
water ice and 23±2% NH 4 SH, where percentages repre¬ 
sent fractional area coverage. This improves the fit in 
the absorbing region to x 2 /N(ABS) = 1.1, a factor of 8 
better than provided by pure NH 3 (see Fig. [12] caption). 
A big improvement is also obtained using the conser¬ 
vative (n=1.4) material in place of NH 4 SH, with a mix 
of 62±3% NH 3 , 35±3% water ice, and 13±2% n=1.4, al¬ 
though the improvement is only by a factor of five in this 
case. (Uncertainties quoted here are deviations that pro¬ 
duce Ay 2 =X%i IN /N , where N is the number of degrees 
of freedom.) The best linear combination without water 
ice produced x 2 /N{ABS) = 5.2, which is a factor of 5 
worse than produced by the best combination with wa¬ 
ter ice. The water ice component reduces the excess I/F 
gradient between 3 /iin and 3.15 /Ain and the excess I/F 
otherwise present between 2.65 /im and 2.85 /iin. The 
third component (either n=1.4 or NH 4 SH) serves to par¬ 
tially fill in the deep minimum produced by pure NH 3 at 
2.96 /im. The linear combination analysis clearly favors 
NH 3 , water ice, and NH 4 SH in that order. Our next step 
was to determine whether what we think is a somewhat 
more realistic physical configuration, i.e. a horizontally 
homogeneous cloud structure, could achieve similar fit 
quality. The red curve in Fig. [T2B is for the last of three 
such structures we tried. 


6.4. Models with composite particles 

Numerical experiments with thick layers of coated par¬ 
ticles containing substantial fractions of NH 3 as the main 
tropospheric haze layer produced excessively strong NH 3 
signatures unless the fraction of NH 3 was so small that 
the desirable weaker features were also eliminated. We 
found a more successful strategy was to split the main 
tropospheric haze layer into two parts, one basically con¬ 
servative, and one dominated by NH 3 . There are two rel¬ 
atively plausible implementations of this approach that 
are both capable of suppressing the stronger NH 3 fea¬ 
tures while retaining enough of the NH 3 character to 
match the VIMS observations. The first (Model 1) is to 
put an optically thin layer of NH 3 particles above a layer 
with relatively flat spectral properties in the 3-/Am region. 
The second (Model 2) is to put an optically thin layer of 
conservative particles above an optically thick layer of 
particles containing a significant component of NH 3 . In 
the following we consider both options and third (Model 
3) that is a variant of the second. The model spectra are 
compared to the measured storm spectrum in Figs. [l3l[T5l 
with model vertical structures displayed in Fig. [TG] The 
complexity of these models, involving different materials, 
different radii, different core fractions and combinations 
of core and shell materials, as well as optical depth and 
pressure of each layer, is so large that we cannot demon- 
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strate uniqueness for any of these solutions. 

6.4.1. Model 1: an NH 3 layer on top. 

The first model we considered consists of a deep water 
cloud at 1.2 bars, a main haze sublayer (haze2 in Fig. 
[HI legend ) of composite particles extending between 210 
mb and 350 mb, an NH 3 haze layer from 110 mb to 210 
mb, and a stratospheric haze that is of little physical sig¬ 
nificance because of uncertain I/F measurements at low 
I/F values. The parameters of this model are defined in 
Table[4] In Fig.IHlthe corresponding model spectrum is 
compared to the spectrum from the head of the storm. 
In Fig. [HR we see that the model spectrum does a good 
job of capturing the main features of the measured spec¬ 
trum, including the spectral features between 2.9 and 
3.1 /im, which are due primarily to NH 3 ice. No other 
plausible material in Saturn’s atmosphere can produce 
the sharp feature near 2.96 /rrn or the secondary bump 
at 2.05 pm. The NH 3 layer has a relatively small optical 
depth (t=1) so that the underlying bright layer can limit 
the depth of the sharp absorption feature at 2.96 pm. 
These NH 3 ice particles are assumed to be spheres of ra¬ 
dius 1.5 pm. We found that significantly larger particles 
tended to produce to large an I/F at longer wavelengths, 
while smaller particles made it more difficult to match 
spectral features near 3 pm. 

The underlying layer that provides most of the re¬ 
flected light visible outside the atmosphere is here com¬ 
posed of composite particles. After trying fits with spec¬ 
trally flat particles in this layer, we were unable to keep 
the I/F at 3.1 pm low enough without causing problems 
at other wavelengths. An example is shown in Fig. IHB. 
which shows the effects of several different modifications 
to the nominal compositional model. In one case we used 
a homogeneous particle of refractive index n=1.4 (com¬ 
puted as a composite particle in which both core and 
shell are the same), which produces the spectrum shown 
by the dot-dash curve. This is seen to be too bright at 
2.9-3.15 pm, while the fit over the rest of the spectrum is 
very good. The nominal model, which uses a small core 
of H 2 0 (40±5% of the total radius), takes care of this 
problem nicely. We also tried to use NH 4 SH as a shell 
material (the dashed curve in Fig. IHI 1). but this causes 
problems at many places in the spectrum. Although with 
some minor adjustment it could produce a reasonable fit 
in the 2.9-3.1 pm region, it tends to be too bright at 1.9- 
2.0 pm, at 2.85 pm, and beyond 3.9 pm, as well as being 
too dark near 3.7 pm because of an absorption by NH 4 SH 
in this region. Some of these problems can be fixed by 
adjusting other parameters, but we could not find a fix 
for all of them simultaneously. It is much easier to fit the 
spectrum using a shell composed of the n=1.4 material 
that seems to be the main component of clouds in other 
(non-storm) regions of Saturn’s atmosphere. 

We also considered a different composition for the deep 
and optically thick cloud that is mainly needed to block 
thermal radiation beyond 4.8 pm. Using the n=1.4 mate¬ 
rial for that cloud produces too high an I/F in the 2.7-2 .8 
pm region as well as allowing too much thermal radia¬ 
tion to leak through near 5 pm. To fix the thermal leak 
would require either greatly increased optical depth or an 
additional absorption. Water fixes both these problems 
because it has absorption in both regions of the spec¬ 
trum, while still providing the needed high reflectivity 


for the shorter wavelength windows at 1.26 pm and 1.62 
pm. 

6.4.2. Model 2: conservative particles on top. 

Perhaps a somewhat more realistic physical model is 
the alternate structure summarized in Table [5j with the 
corresponding spectrum displayed in Fig.[H] In this case 
the top part of the main cloud layer consists of a pure 
conservative substance, which we assume to be the same 
material seen outside of storm regions and is assumed to 
have a refractive index of n=1.4+0i. The NH 3 appears 
in the lower part of the main cloud layer as a coating on 
a small core of H 2 O (also about 40% of the total radius). 
The deep layer is composed of pure H 2 O ice, as in the pre¬ 
vious composite model. This produces a spectral match 
which is comparable to the first model, but is more easily 
understood physically. If the storm convection begins at 
the water cloud level, then as the water particles rise and 
freeze they will eventually serve as (large) condensation 
nuclei for NH 3 which will eventually reach saturation and 
form a coating on the water cores. This makes more sense 
than forming coatings of n=1.4+0i material, since that 
material is not even known to be formed by condensation 
at the relevant level of the atmosphere. The top layer of 
n=1.4+0i material can perhaps be understood as part of 
the original main cloud layer pushed to higher altitudes 
by the converting tower that brings up water from the 
deep levels of the atmosphere. This layer has the im¬ 
portant spectral effect of filling in the deep minimum at 
2.96 /un that would otherwise be visible due to the layer 
below it, as can be seen from the spectrum in Fig. [HR 
that is obtained when the layer is removed. Although 
we do not seem to need NH 4 SH or find it useful either 
in the form of a pure particle or as a coating or core in 
a two-component composite particle, it is expected that 
particles are likely to contain a water core, a layer of 
NH 4 SH , followed by a layer of NH 3 . Unfortunately, we 
do not have tools for modeling such particles. 

6.4.3. Model 3: conservative top layer with an extended deep 

cloud. 

It is also possible and probably more physically realis¬ 
tic to model the deep cloud as vertically extended instead 
of as a compact layer. Such a model, with a uniform op¬ 
tical depth per bar from the 10 bar level to the bottom 
of the main (haze 2 ) cloud, can indeed fit the observa¬ 
tions just as well as the prior compact model for the deep 
cloud, as long as the cloud density is sufficient to block 
out thermal radiation, which requires about 10-20 opti¬ 
cal depths per bar. We also found that larger particles 
in the NFR-containing layer provided a somewhat better 
match to the spectral features near 3/tm, although main¬ 
taining a good overall fit to the spectrum in this case re¬ 
quired adding some absorption at long wavelengths. To 
accomplish that, we made the previously conservative 
substance absorbing between 4.3 and 5.3 /mr by linearly 
increasing its imaginary index from zero to 0.03 over that 
region. We also decreased the particle radius in the up¬ 
per sublayer from 1.5 /tm to 1.0 /tm to further reduce its 
reflectivity at longer wavelengths. The model parame¬ 
ters are provided in Table [ 6 ] Some of these were fixed 
manually, while others with uncertainty estimates were 
adjusted to minimize x 2 within the assumed constraints. 
This revised model is compared with observations in Fig. 
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TABLE 4 

Fit parameters for composite storm cloud structure Model 1. 


Particle 

Layer boundaries 

Particle 


Composition 

Ptop 

Pbottom 

radius 

Optical depth 

nh 3 

110 mb 

210 mb 

1.50 ^m 

1.0 

n=1.4+0i shell 
on H 2 O core 

210 mb 

350 mb 

1.5 fim 

with 

0.6 fim 

8.0 




core 


h 2 o 

1000 mb 

1200 mb 

2.0 fim 

80 


We also included a stratospheric haze between 0.5 mb and 2 mb, which here was treated as a 
double-HG scatterer with t=0.01 at 2 pm, and a A -3,5 wavelength dependence. Nominal gas 
mixing ratios and profiles were assumed as described elsewhere, except that the PH 3 scale 
height was assumed to be 50% of the gas scale height. 



Fig. 13.— A: VIMS measured spectrum of the storm head (thick gray line) in comparison with a complete model spectrum (solid line), 
and model spectra with different layers removed one at a time (see text for discussion). B: Model spectrum (solid curve) compared to 
spectra for modified models in which the deep water cloud is replaced by a cloud with n=1.4 (dotted), or the H 2 O core of the 210-350 
mb haze is replaced by an n=1.4 core (dot-dashed), or the n=1.4 shell of the 210-350 mb haze particles is replaced by an NH 4 SH shell 
(dashed). See text for discussion. 


Q! and also with the best linear combination heteroge¬ 
neous model in Fig. [T3] (red curve). While this provides 
better fit to the spectrum near 3 /im than the other hor¬ 
izontally homogeneous models, it is still does not quite 
match the quality of the best heterogeneous model. This 
probably means that there is more physics needed in the 
homogeneous model to capture more complexity in com¬ 
position an d vert ical structure. 

iSanz-Reciuena et ahl (l2012h analyzed ground-based vi¬ 
sual imaging observations covering UV and a range of 
methane bands for the purpose of characterizing the 
storm’s cloud structure. Their observations were mod¬ 
eled with a 3-layer structure consisting of a stratospheric 
haze, a tropospheric haze, and a bottom cloud. They 


found that the tropospheric haze optical depth was lower 
in the storm cloud than it was in the non-storm region 
(15±5 vs 25±5), although the opacity was concentrated 
at higher altitudes as they inferred a cloud boundary 
change to 100 - 300 mb (storm) from 100 to >1000 mb (no 
storm). This concentration of optical depth at low pres¬ 
sures is at least roughly consistent with our inferences 
from near-IR observations, and also with our finding of 
a more extended main layer in the non-storm regions. 
However, as shown in Fig. 1161 we find that the storm re¬ 
gion is better fit with a two-component main layer, with 
most of the optical depth in the lower sub-layer, and 
with the upper sub layer extending to somewhat lower 
pressures than found in the non-storm regions. 
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Fig. 14.— A: VIMS measured spectrum of the storm head (thick gray line) in comparison with a complete model spectrum (solid line) 
for the alternate composite model defined in Table [5] model spectra with different layers removed one at a time are als o sh own (see text 
for discussion). The removal of the stratospheric haze layer (not shown) has essentially the same effect as shown in Fig. 1131 B: Alternate 
model spectrum (solid curve) compared to spectra for modifications in which the deep water cloud is removed (dotted), or the H 2 O core 
of the 230-300 mb haze is replaced by an NH 4 SH core (triple dot-dashed). See text for discussion. 


TABLE 5 

Fit parameters for storm cloud composite structure 
Model 2. 


Particle 

Composition 

Layer boundaries 
Ptop Pbottom 

Particle 

radius 

Optical 

depth 

n=1.4+0i 

150 mb 

230 mb 

1.50 fim 

1.1 

NH 3 shell on 
H 2 O core 

230 mb 

300 mb ( 

1.5 fim with 
0.6 fim core 

9.0 

h 2 o 

900 mb 

1000 mb 

2.0 fim 

80 

Stratospheric haze and PH 3 scale height 

as in Table [4] 


TABLE 6 

Fit parameters for storm head Model 3, which includes a 

VERTICALLY EXTENDED WATER CLOUD. 

Particle 

Composition 

Layer boundaries 
Ptop, mb Pbot, mb 

Particle 

radius 

Optical 

depth 

n=1.4+0i 

140±79 

189±25 

1.0 fim 

1.46F0.06 

NH 3 shell on 
H 2 O (core) 

189±25 

544±135 

2.25 /mi 
(0.84 fim) 

7.1F1.3 

H 2 O 

544±135 

10 bars 

2.0 fim 

20 /bar 


Stratospheric haze and PH 3 scale height as in Table [4] 


7. SUMMARY AND CONCLUSIONS 

From an analysis of VIMS spectra of the giant Saturn 
storm of 2010-2011 and its surroundings, we have reached 
the following conclusions. 

1. VIMS spectra outside the storm region exhibit a 
high degree of uniformity in the fractional depth of 
a 2.95-3.0 (im absorption feature across the Saturn 
disk, with essentially no dependence on latitude or 
view angle, implying that the feature is primar¬ 
ily due to a responsivity error in the vicinity of a 
2.98-/xm joint between order-sorting filters. It also 
has. a slight variation with the x coordinate of the 
VIMS images, but essentially no variation with the 
y coordinate. We developed an accurate correction 
algorithm that removes this feature, allowing us to 
detect true absorption by NH 3 in the vicinity of 
storm clouds. 

2. Most near-IR spectra of Saturn outside the storm 
regions can be very well fit with relatively simple 
models consisting of ( 1 ) a deep and often optically 
thick cloud layer that serves to block thermal ra¬ 
diation in most places, but is thin enough, or deep 
enough, in some regions to allow 5-/im detection of 
atmospheric composition and horizontal structure 
at the several bar level, ( 2 ) a middle “main” cloud 
layer extending from a top pressure of 111-178 mb 
to a bottom pressure of 577-844 mb, with optical 
depths of 3.5-8 .6 at 2 /jm, and composed of parti¬ 
cles that are relatively conservative scatterers, with 
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Fig. 15.— A: VIMS measured spectrum of the storm head (thick gray line) in comparison with the model 3 spectrum (solid line), and 
modified spectra with different values for the optical depth per bar of the vertically extended deep water cloud. B: Model 3 spectrum (solid 
curve) compared to spectra for modified models with slightly different parameter choices made one at a time. See text for discussion. 



Fig. 16.— Comparison of non-storm cloud structure models A and C (from Fig. [8]) with Models 1-3 for the head of the storm. The two 
X 2 values given for Models 1-3 are for the 1.268-5.15 (i m region (first) and the 2.65-3.2 fun region (second), with numbers of comparison 
points being 251 and 51 respectively. 


an acceptable refractive index of n=1.4+0i, and (3) 
a stratospheric haze of low optical depth, but which 
cannot be accurately constrained by near-IR VIMS 
observations because they are too uncertain (and 
generally too large) at very low I/F values (near 
the 1 DN level). 


3. Outside the storm-affected regions, there is no ev¬ 
idence for particulate absorption at 3 fi m, further 
confirming the peculiar result that at least the up¬ 
per several optical depths of Saturn’s main visible 
cloud layer is not made of NH 3 , NH 4 SH, or N 2 H 4 . 
We cannot currently rule out P 2 H 4 , but its pres- 
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ence seems likely to be testable in the future from 
measurements near 4.3 fim once optical constants 
for solid diphosphine are determined. 

4. Inside the storm head we found evidence for a 
prominent 3-/mr absorber in Saturn’s storm clouds 
but little indication of an absorber active near 2 
/jrn. A weak 2-/im feature is not too surprising if 
the NH 3 ice is confined to a cloud of relatively small 
optical depth overlying a relatively bright cloud 
that fills in the lack of reflectivity in the absorbing 
regions of the NH 3 ice spectrum, or if the ammonia- 
containing cloud layer is underneath a conservative 
scattering layer that serves the same purpose. An 
uncovered optically thick layer of pure NH 3 ice par¬ 
ticles would produce a noticeable feature at 2 /iin 
as well as a very prominent feature 2.96 /irn, and 
thus a thick cloud of pure NH 3 ice particles is not 
plausible as the dominant upper tropospheric cloud 
in the most absorbing storm regions. 

5. We considered four plausible candidate materials 
that might contribute significant absorption in the 
3-/im region: NH 3 , N 2 H 4 , NH 4 SH, and H 2 O. Fits 
to the observed spectrum outside of the 2.8-3.1 /jm 
region of the spectrum were obtained for a model in 
which each material was taken to be the only com¬ 
ponent of the main cloud layer. We then compared 
the features inside 2.8-3.1 pm part of the spectrum 
as well as outside with the VIMS measured spec¬ 
trum. The only substance that came close to fitting 
the observed spectrum was NH 3 , but it had a very 
sharp absorption feature at 2.96 pm that was not 
present in the observed spectrum. 

6 . We also considered linear combinations of the spec¬ 
tra obtained from the pure composition fits, as 
might be obtained from a heterogeneous horizon¬ 
tal mix of different particle types. In this case we 
found that the 2.96-pm feature could be reduced 
in amplitude by adding contributions from other 
materials. The best fits were obtained with a mix 
of NH 3 , NH 4 SH, and H 2 O, using areal fractions of 
0.55, 0.23, and 0.22 respectively. 

7. The best fits to VIMS spectra using a horizon¬ 
tally homogeneous vertical structure model con¬ 
tained a deep water cloud, perhaps beginning at 
10 bars, but certainly containing significant opacity 
above the 1-2 bar range, a main cloud layer con¬ 
sisting of two sub-layers, one of low optical depth 
(t S3 1 ) above a lower sub layer of significant op¬ 
tical depth (r ss 9) and an optically thin strato¬ 
spheric haze that we could not reliably constrain 


Acton, C. H., 1996. Ancillary data services of NASA’s Navigation 
and Ancillary Information Facility. Planet, and Space Sci. 44, 
65-70. 


because of VIMS limitations at very low signal lev¬ 
els. We initially found two options for the main 
layer produced good fits: (1) particles of NH 3 ice 
in the top layer and particles of n=1.4 material 
(unknown composition also present almost every¬ 
where else on Saturn) coating a core of water ice 
in the lower sublayer, or (2) particles of n=1.4 in 
the top sublayer and NH 3 in the lower sublayer as 
a coating on H 2 O cores, with a core radius that is 
40% of the total radius (1.5 n m), and a core vol¬ 
ume of 6 %. A slightly better fit was obtained with 
option ( 2 ) using smaller particles with some long¬ 
wave absorption in the upper main cloud layer and 
somewhat larger NH 3 -coated water ice particles in 
the lower main sub layer. 

8 . Although NH 4 SH was the best third component of 
our linear combination spectra that represent hori¬ 
zontally heterogeneous mixtures, we did not find a 
useful role for it in our horizontally homogeneous 
composite models. It might be present in small 
amounts in several layers, but when it becomes a 
significant fraction of the material in the particles, 
it presents spectral problems that are much easier 
to solve with other materials. However, it is pos¬ 
sible that NH 4 SH might prove useful as part of a 
three-component particle (core with two coatings), 
which we do not have the tools to model. 

9. Both heterogeneous and homogeneous models of 
the VIMS storm spectra strongly favor contribu¬ 
tions by water ice, which is the first strong spec¬ 
troscopic evidence for water ice in Saturn’s at¬ 
mosphere, found near the level of Saturn’s visible 
cloud deck where it could only be delivered by pow¬ 
erful convection originating from ~200 km deeper 
in the atmosphere, as expected from numerical dy¬ 
namical models. 

10. The increased scale height we found for PH 3 in 
the storm head relative to the surrounding regions 
(0.5 compared to 0.375 times the pressure scale 
height) is consistent with increased upward con¬ 
vective transport inside the storm. 

11. We find strong evidence for the presence of NH 3 
ice at the top of the clouds, but at least some 
(and perhaps much) of the cloud seems to be of 
the same relatively conservative (at near-IR wave¬ 
lengths) material that forms the main component 
of clouds over the non-storm regions. 
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